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ABSTRACT 

We present a study of the morphological nature of redshift z ~ 0.7 star-forming galaxies using 
a combination of HST/ACS, GALEX and ground-based images of the COSMOS field. Our sample 
consists of 8,146 galaxies, 5,777 of which are detected in the GALEX near-ultraviolet band (2310A or ~ 
1360A rest-frame) down to a limiting magnitude of 25.5 (AB), and all of which have a brightness of 
F81AW (HST) < 23 mag and photometric redshifts in the range 0.55 < z < 0.8. We make use of the 
UV to estimate star formation rates, correcting for the effect of dust using the UV-slope, and of the 
ground-based mutli-band data to calculate masses. For all galaxies in our sample, we compute, from 
the ACS F814W images, their concentration (C), asymmetry (A) and dumpiness (S) as well as their 
Gini coefficient (G) and the second moment of the brightest 20% of their light (M20). We observe a 
bimodality in the galaxy population in asymmetry and in dumpiness, though the separation is most 
evident when either of those parameters is combined with a concentration-like parameter (C, G or 
M20). We further show that this morphological bimodality has a strong correspondence with the 
FUV - g color bimodality, implying that UV-optical color predominantly evolves concurrently with 
morphology. We observe many of the most star-forming galaxies to have morphologies approaching 
that of early-type galaxies, and interpret this as evidence that strong starburst events are linked to 
bulge growth and constitute a process through which galaxies can be brought from the blue to the 
red sequence while simultaneously modifying their morphology accordingly. We conclude that the 
red sequence has continued growing at z < 0.7. We also observe z ~ 0.7 galaxies to have physical 
properties similar to that of local galaxies, except for higher star formation rates. Whence we infer 
that the dimming of star-forming galaxies is responsible for most of the evolution in star formation 
rate density since that redshift, although our data are also consistent with a mild number evolution. 
Subject headings: galaxies: evolution — galaxies: fundamental parameters (masses, morphologies, 
radii, star formation rates) — surveys 
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1. INTRODUCTION 

Star formation has been in decline, in the Universe, 
for the past 8 billion years. This disc overy, manifested 
in the now well-known Madau diagram ULillv et al.ll 19961 
iMadau et alJll996[ ). has been a remarkable culmination 
of last decade's research in galaxy evolution. The details 
of this decline remain, however, surprisingly elusive: the 
reason lying in the complex nature of star formation it- 
self. Fuel exhaustion (in situ gas consumption), reduc- 
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tion in merger rate and environmental effects could all 
share the responsibility. As many physical character- 
istics of galaxie s correlate with their Hubble type (e.g. 
lKennicu"ttlll998f) . their classification within redshift sur- 
veys has been a natural segue in the investigation of the 
decline of star formation. 

Color, spectral class or morphology are among the 
most commonly used criteria to separate galaxies and 
were quick l y app lied by investigators. For example, 
iLillv et al] (|1995f) found that the luminosity density 
of blue galaxies brightens substantially up to redshift 
z ~ 1, whereas that of red galaxies does not, a result 
confi rm ed by toda y 's mu ch larger samples (jFaber et al.1 
I2005H . lEllis et alj (|1996l) divided galaxies in the Aut- 
ofib/LDSS sample according to their [Oil] equivalent 
width and found a strong evolution in the volume den- 
sity of moderate to low-luminosity o bjects with strong 
[Oil] emission while iHevl et ail (|1997l ) divided the same 
sample into spectral types and found that late-type spi- 
rals were the ones dominating the evolution of the blue 
luminosity function. 

Attempts at resolving distant star-forming galaxies 
into morphological types date back to the problem of 
faint blue galaxies (see lEllisI 119971 for a review) for 
which IBrinchmann et all ( 1998t ) concl uded, using HST 
pointings a t locations in th e CFRS ([Lilly et al.l Il995h 
and LDSS (jEllis et al.l ll996f) fields, that peculiar galax- 
ies, identified to be r esponsible for the faint blue excess 
(jGriffiths et al.lfl99l . were also the main cause for the 
rapid evolution of the blue luminosity functi on observed 
in tho se redshift surveys. Recently, however. IWolf et al.l 
(2005), with a sam ple of 1483 gala xies at z ~ 0.7 ex- 
tracted from GEMS (|Rix et alJl2004f ). found spiral galax- 
ies to actually dominate the overall UV(2800A) lumi- 
nosity at that redshift (though with irregular galaxies 
still being prevalent at faint magnitudes), implying that 
their fading, accompanied with a similar migration to 
lower UV-luminosities in irregulars, must lead the de- 
cline of star formation density, rather than a decrease in 
merging rate. Spitzer 24fim observations of that same 
sample show that most of the IR-emission associated 
with dust-reprocessed UV-light from young stars, emis- 
sion that declines with redshift even fast er than that of 
the escaping UV (iLe Floc'h et al.l 120051). is also domi- 
nated by spiral galaxi es (Bell et al.M2005[ ). Moreover, 
iMelbourne et al.l (|2005f ) compared LIRG morphologies at 
rcdshifts 0.1 < z < 1 and confirmed that high-redshift 
(z > 0.5) LIRGs are dominated by spirals unlike low- 
redshift ones (z < 0.5) which are mostly peculiars. They 
interpret that as a depletion of gas supply causing spirals 
to fall to sub-LIRG levels of star formation, while pe- 
culiar morphologies, characteristic of mergers, continue 
to experience strong bursts of s tar fo rmation. Along 
the same lines. iMenanteau et al.l (|2006h . using the paral- 
lel NICMOS observations of the UDF, also showed spi- 
rals to dominate Psfr at all redshifts up to z > 1, at 
which point irregular/peculiar galaxies, which show the 
sharpest rise, become equally important. These many 
evidences suggest that spirals play a crucial role in the 
last 8 Gyr evolution of galaxies. 

Meanwhile, a wealth of information about low-redshift 
galaxies has al so emerged. In p articular, large surveys 
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gaiaxies nas al so emerged, in p ; 
such as SDSS (|York et all [20001 



shift Survey (e.g. iMadgwick et al.ll2002D have really pro- 
pelled statistical studies of galaxy properties. Although 
trends of c olor vs. morphology h ave been know n for a 
long time (|de Vaucouleurslll961f ). IStrateva et~ail 1)20011 ) 
showed that the color distribution of galaxies was not 
smooth, but doubly-peaked with early-types dominat- 
ing the red population and late-types the blue one. 
iKauffmann et all (|2003af) also observed this bimodality 
in the £>„(4000) index, indicative of a division in the 
population between galaxies dominated by an old stellar 
population and ones that experienced re cent episodes of 
star formation. Kauffmann et al. (2003b) further showed 
that this separation occurs at M* ~ 3 x 1O 1O M0 and that 
lower-mass galaxies besides having young stellar popula- 
tions also have disk- like structural parameters, whereas 
higher-mass ones have old stellar populations a nd bu lge- 
like structural parameters. IBrinchmann et ail (|2004l ) re- 
inforced these conclusions by directly calculating specific 
star formation rates (SFR per unit mass) for ~ 150, 000 
galaxies, and observed the same divide of high-sSFR 
and low-sSFR wit h age, mass and struc tural parame- 
ters/Hubble type. iBlanton et alj (|2003bl ) labeled these 
two pop ulations the red and the blue sequences. More 
recently, iBell et ail (|2004h found the same color bimodal- 
ity at a ll redshifts up to z ~ 1 using the COMBO-17 
sample (|Wolf et al.ll2003fl. 

With COSMOS (|Scoville et al.ll2007bft . it is now possi- 
ble to expand statistical analyses of higher redshift galax- 
ies to levels that allow comparison with SDSS or 2dF. In 
this paper, we combine HST/ACS and GALEX coverage 
of th e COSMOS field (jKoekemoeri [2007t ISchiminovichl 
l2007f ) to study the morphological properties of star- 
forming galaxies at z ~ 0.7. In view of the streneous task 
of classifying large numbers of objects by eye and of the 
subjectivity it carries, for this study as well as for future 
comparison we chose to follow the pa th of automated 
classification. lAbraham et al.l (1996ba) demonstrated 
the usefulness of such an approach by using concentra- 
tion and asymmetry measurements to classify galaxies 
in the MDS and HDF ac cording to their l ocatio n in the 
C-A plane. In parallel, Odcw ahn et al.l (|1996h used a 
neural-network code on the HDF sample and arrived 
at results similar to lAbraham et ail (|1996af ). Following 
along this path, people have devised and applied several 
ingenious algorithms for morphologi cal classificat ion, the 
most popular being the Sersic index (persid [19^ bulge- 
to-disk decomposition iSiinanl '199 a; iPeng et all 12002) . 
and shapelet decomposition (Refregier 2003). These, 
however, all require fitting objects to a set of parame- 
ters or functions. In this paper, we pre ferred to follow 
the pu rely mensurational approach of lAbraham et al.l 
(1996b]), an d expanding it by adding measurements of 
dumpiness (|Conselicej [2003) as well a s of the recently 
developed Gini and M20 coefficients (jLotz et al.l 120041 : 
lAbraham et al.l l2003f ) to use as our morphological pa- 
rameters, and classification criteria. 

In this paper, we present morphological characteristics 
of the galaxy population at z ~ 0.7, and study their 
relation to physical parameters. We also compare both 
their morphological and physical properties with those of 
low-redshift galaxy samples in the literature. We focus 
most of our attention on the relation between star for- 
mation rate and morphology, and interpret our results in 
the framework of galaxy evolution. We discuss implica- 
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tions for evolution scenarios since z ~ 0.7 in the context 
of the literature, with an emphasis on blue to red se- 
quence evolution. The paper is organized as such: we 
first shortly describe, in section 2, the COSMOS survey, 
its data, as well as the GALEX observations and data, 
before discussing our sample selection. We outline our 
approach to morphological analysis in section 3, present 
the results of our investigation in section 4, and discuss 
their interpretation and implications in section 5. 

2. OBSERVATIONS AND DATA 
2.1. EST/ ACS Observations 

We make use of the full HST coverage of the COS- 
MOS field, which consists of 542 HST/ACS images 
with depth of I < 27 mag (AB,10cr point source), 0.09" 
FWHM resolution (with 0.05" pixels) and whose mosaic 
spans an area of 2 deg 2 . An overview of the COSMOS 
project is given in lScoville et al.1 (|2007b| ) with detail de- 
sc ription of the ACS o bserv ations and data r eduction 
in IScoville et~aT1 (|2007ah and iKoekemoerl (|2007h respec- 
tively. 

2.2. Ground-based Observations and Catalog 

Ground-based follow-up observations have been per- 
formed using the CFHT (u* and i bands), Sub- 
aru/SuprimeCam (BVgriz), Kitt Peak/Flamingos (K- 
band) and CTIO (also K-band) telescopes, providing 
deep coverage, with typical limiting magnitudes of 27 
(AB, 3tr), of the field from the u to z bands (mf imit = 
25.8), as well as shallower imaging in the K-band 
( m Hmit = 21.6). Details of the ground-b ased observa- 
tions and data reduction a r e pres ented in ICapak et al.l 
(|2007ft and iTaniguchi et alj (120071) A m ulti-wavelength 
photometric catal og dCapak et all 120071) w as generated 
using SExtractor (|Bertin &: Arnoutall996h . with the i- 
band as the selection wavelength. We further performed 
SED fitting of this multi-band data and calculated pho- 
tometric redshifts for galaxies with i < 25 mag (AB) 
(jMobasher et al.ll2007f ). Our photometric redshifts have 
an rms of a((z phot - z spec )/(l + z spec )) = 0.031 with 2% 
outliers. We make use of these photometric redshifts in 
selecting our sample. 

2.3. GALEX Observations 

We used GALEX ([Martin et al.ll2005al ). which has a 
circular field-of-view of 1.2° in diameter, to observe the 
COSMOS region in ultraviolet light with four pointings 
of ~ 50 ks each. These observations, performed as part 
of the GALEX Deep Imaging Survey, reach a limiting 
magnitude of ~ 25.5 mag (AB) in the near-ultraviolet 
band (NUV). GALEX has a resolution of 5.6" in the 
NUV which corresponds to ~ 40 kpc at z ~ 0.7, larger 
than the typical size of gala xies. The full details of the 
observations can be found in [Schiminovich (2007). 

Since standard pipeline processing of deep fields can 
sometimes blend two objects into a single detection, 
we employed a different method of source extraction. 
As the vast majority of our sources appear unresolved 
to G ALEX, we de cided to use the DAOPHOT soft- 
ware (Stetson 1992) to measure photometry. We further 
made use of the ground-based COSMOS catalog to feed 
DAOPHOT with position priors. Because of small as- 
trometric offsets, our first step was to align the priors 



with each GALEX image. We ran the phot routine (a 
routine that performs aperture photometry) with cen- 
troid recentering of the objects. We then fit the center 
of the distribution of shifts in the x and y-directions as 
well as in the angle of rotation 9 around the center of 
the pointing, and applied the mean shift to all positions 
obtained from the astrometry. This small (< 1 pix) con- 
stant uniform shift does a good job at realigning position 
priors with objects, and is therefore the only astrometric 
correction we applied. 

We then followed the standard DAOPHOT procedures 
of running phot (performing aperture photometry), psf 
(modeling the psf) and allstar (performing psf-fitting 
photometry) to obtain UV-fluxes for our objects. This 
time we did not allow for recentering in the phot pro- 
cedure, but we did in allstar, since, by looking at the 
residuals, we found psf-fitting to be much better when 
recentering was allowed. The drawback, however, is that 
priors located in regions with no apparent UV were often 
moved, in the process of recentering, to fit a neighbor- 
ing object. The measurements were therefore rematched 
to objects in the original catalog that were located clos- 
est, but no farther than 3", to the measured GALEX 
positions. Lastly, we created masks to eliminate vari- 
ous artifacts, as well as a few very bright stars, in our 
four NUV images, and nulled all detections found inside 
masked regions. 

2.4. Sample Selection 

We aimed at extracting a sample in a narrow range of 
redshifts around z ~ 0.7 bright enough to study mor- 
phology in the ACS images. The choice of redshift 0.7 is 
convenient in that the observed NUV-band roughly cor- 
responds to the (z = 0.1)-frame FUV. This minimizes 
K-corrections and allows for easy comparison with lo- 
cal samples such as SDSS. A narrow redshift range fur- 
ther allows us to obviate the need for morphological K- 
corrections. The selection procedure we employed is the 
following: we first ran our morphological analysis ex- 
clusively on objects with / mag < 23 as morphological 
parameters become less reliable for fainter objects. We 
then removed all objects with (petrosian) radii smaller 
than 0.2", since they are below our resolution limit. As 
those objects are mostly stars and QSO's (figure[T]), this 
cut does not introduce any bias in our sample. Our mag- 
nitude cut does, however, progressively bias us towards 
higher surface brightness objects as we move to smaller 
radii. This effect shows up later in some of our analysis, 
and is discussed in context. We further cleaned our sam- 
ple of stars (and stellar-like object) by applying a cut in 
the r pe tro — I Acs plane (these parameters are described 
in the next section), where stars and galaxies segregate 
unmistakably (figure [T|), before proceeding to remove ob- 
jects that, after visual inspection, were found to be false 
detections or that showed various problems with their 
segmentation (we describe our segmentation technique in 
section 3.1). These fews steps not only clean our sample 
of stars, but also, at the same time, of most but the weak- 
est AGNs. We lastly selected for this study only objects 
with photometric redshifts in the range 0.55 < z < 0.8. 
This redshift bin width translates into a difference of lu- 
minosity of 0.8 magnitudes for objects with the same 
brightness located on both ends of the redshift range. 
In the end, our sample contains 8, 146 galaxies, 5, 777 of 
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Fig. 1. — Petrosian radius vs /^cs-band magnitude for all 
objects in our morphological sample. The dashed-line represents 
our star /galaxy separation, stars lying below the line and galaxies 
above. 



which are detected in the NUV with GALEX. We thus 
detect, in the UV, about 70% of objects with / mag < 23 
and redshift z ~ 0.7. Throughout this paper, we some- 
times utilize, where appropriate, our UV-detected sample 
only, but otherwise normally refer to our full sample of 
8, 146 galaxies. 

2.5. Star Formation Rates and Masses 

We performed our own SED analysis on the combined 
GALEX + ground-bas ed photometric data using the 
KCORRECT software (jBlanton et alJ l2003aT >. and ex- 
tracted from it K-correction estimates for each galaxy. 
We chose such an approach because KCORRECT is de- 
signed to extract the most physically realizable SED by 
using linear combinations of four spectra that are charac- 
teristic of physical states of galaxies, from intensely star- 
bursting to quiescent. We applied the K-corrections to 
u* and NUV-band photometry to obtain restframe FUV 
and NUV absolute magnitudes from which we derived a 
UV-slope which, given the relation betwee n (the UV- 
slope) and Apuv (the FUV attenuation) (jSeibert et alJ 
2005), provides us with a dust correction factor. 17 We 
then converted the corre cted FUV-lumino sities to star 
formation rates using the iKennicuttl (j!998l ) relation be- 
tween star formation rate and UV-continuum luminosity: 

SFR(M Q year" 1 ) = 1.4 x 10" 28 L„ uv (ergs s" 1 Hz" 1 ) 

(1) 

For objects without NUV counterparts, restframe FUV 
magnitudes have been derived directly from the fitted 



17 Some observations (Seib ert et al.i 2005; C ortese et aT]|2006l) 
suggest that the Apjjv — & relation is steeper in starburst galaxies 
than in normal galaxies. Since we are using the relation for nor- 
mal galaxies, it is possible that, for star bursts, our star formation 
rates are slightly underestimated. This, however, would not affect 
the qualitative behavior of star formation rate in relation to other 
properties, which is what we focus on in this paper. 



SEDs. Furthermore, because of the degeneracy between 
age and dust in red galaxies, and because of the uncer- 
tainty in FUV magnitudes derived for objects with no 
GALEX detection, we decided not to use (3 as a proxy for 
dust attenuation for these galaxies, but instead, to sim- 
ply apply a moderate constant dust correction of 4-0.5 
in log SFR (which is equivalent to a Apuv 01 1-25 or 
an E(B-V) of 0.151). This is a reasonable correction for 
early-type galaxies which constitute most of these ob- 
jects. On the other hand, this method completely misses 
the most heavily obscured galaxies, such as could be some 
ULIRGs, though these are far less common. The up- 
coming Spitzer data release for the COSMOS field will 
be extremely helpful in the study of these objects. For 
now, we need to leave those with UV-fluxes below our 
detection limit behind, i.e. with star formation rate es- 
timates in the range of quiescent galaxies, much below 
their true value. Because of the discrepancy in the qual- 
ity of our measurements between UV-detected and non- 
UV-detected objects discussed above, we clearly differ- 
entiate the two populations in our plots and analysis. 

After conversion to a star formation rate, our limiting 
magnitude of tunuv = 25.5 corresponds to log SFR = 
0.11 or a star formation rate of about 1 M Q yr _1 for a 
z = 0.7 galaxy with an Apuy of 1.25. Throughout this 
paper, we thus also refer to objects in our UV-detected 
sample as star-forming galaxies. 

We also applied our derived K-corrections to obtain 
res tframe B and U - band absolute magnitudes and used 
the lBell fc de Jond (|2001h relation between B - V color 
and the ratio of mass to U-band luminosity 

log(M*/Ly) = -0.734 + 1.404 x (B - V) (2) 

with a scaled Salpeter IMF to calculate masses for 
all objects in our sa mple. The scale Salpeter IMF 
(|Bell fc de Jond |2001| ) has a shallow er slope at low 
masses, similar to the Chabrier IMF (|Chabrierl l2003f) . 
This method is accurate to about 0.1-0.2 dex with most 
uncertainties coming fro m bursts of star forma tion, dust 
and uncertainties in the lBell fc de Jonel (|2001l ) models. 

3. MORPHOLOGICAL ANALYSIS 

Given the large nature of our sample, it is impor- 
tant to use both an automated and consistent morpho- 
logical classific ation scheme. We there f ore chose to fol- 
low th e work of | Abraham et al.l (Il996bf). IConselice et al.l 
(120001). IConselicd (|2003ft . lAbraham et all (|2003f ) and 
lLotz et all (|2004f ) and use their non-parametric ap- 
proaches, thus computing Concentrati on(C), Asymme - 
try(A) and Clumpiness(S) parameters (|Conselicd l2003). 
as well as the Gini coefficient (G) and second order mo- 
ment of th e distribution of the brightest 20% of the 
light(M20) (|Lotz et al.ll2004h for all objects in our sam- 
ple. We briefly describe these parameters in the following 
sections, but refer the reader to the papers cited above 
for a full description. 

3.1. Size and Segmentation 

Before measuring any morphological parameter, one 
needs to assign a reg ion of the image to every galaxy. The 
standard approach (Conselicc 2003) has been to take a 
circular aperture of radius 1.5 x r petro , where r petro , the 
Petrosian radius, is the radius at which rj{r) = 0.2, and 
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where rj(r) is defined as 
7?(r) 



(3) 



that is the ratio of the surface brightness at a given ra- 
dius to the mean surface brightness within that radius. 
Over a fixed surface brightness cut, this method has the 
advantage of being far less affected by surface brightness 
dimming. 

The Gini coefficient and the second order moment of 
the light, h owever, require a full segmentation of ev- 
ery object. lLotz et all (|2004f ) use the isophote of sur- 



face brightness fi = /u(r. 



petro / 



on a smoothed version 



of the image, the smoothing kernel being a gaussian 
of cr = 0.2 x r petro , as the boundary of their galax- 
ies. We followed the same prescription, though with a 
tophat smoothing kernel of diameter 0.3 x r petro . We 
also imposed a minimum surface brightness of [i m in 
0.6 x aback ground (after background subtraction). By 
summing all the flux within our segmentation maps, we 
were able to estimate the total -Oics-band flux for each 
galaxy, which we then used to normali ze all of our mor - 
phological and size parameters. Unlike IConselicel (2003) 
who used thumbnail images, our algorithm extracts ob- 
jects from larger images, and we therefore decided to also 
extend the use of segmentation maps to the computation 
of asymmetry and dumpiness since, compared to circu- 
lar apertures, they are less likely to pick up light from 
neighboring objects. 



Apart from 



petro 



we also make use of three other size 



measurements: r 2 o%, r 50 % and rg %. These represent 
the radii encompassing respectively 20, 50 and 80% of 
the total /^cs-flux of the galaxy, and are obtained by 
summing the flux inside a circular aperture of expand- 
ing radius until the respective percentages of the light 
are attained. We utilize r 2 o% and r 80 % in the computa- 
tion of concentration, and we use r 50 % (converted to a 
physical scale) as our parameter for size. Because our 
value of r 50 % is measured with respect to the total flux 
inside the Petrosian radius, we systematically underes- 
timate its actual value. Therefore, in order to attempt 
to recover to true size of our objects, we multiplied our 
values of r 50 % by a correction factor. This correction fac- 
tor was calculated by, first co mparing tru e values of r 50 % 
for theoretical Sersic profiles (|Sersid 1968) with bcrsic in- 
dices of 1,2,3 and 4, to the radius at 50% of the flux inside 
their Petrosian radius, then calculating the concentration 
values for those four profiles (there exists a one-to-one 
monotonic relation between Sersic index and concentra- 
tion), and finally obtaining a fit for 
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as 



a function of concentration. The result is shown in fig- 
ure [21 The points represent the theoretical calculations 
for Sersic profiles of Sersic indices 1,2,3 and 4, and the 
line represents the fit to the relation. We also limit the 
correction factor to 2.0, hence all objects with C > 4.65 
simply have their r 50 <% doubled. 

3.2. Concentration 

Concentration is defined by the ratio of the radius con- 
taining 80% (rgo%) to the radius containing 20% (r 2 o%) 
of the total light. 

r 80% 



C = 5 x loe 



?"20% 



(4) 
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Fig. 2. — Actual 50%-radius over the measured 50%-radius as of 
function of measured concentration, for theoretical Sersic profiles 
of integer indices 1 to 4. The line is a fit through the four calculated 
points and represents the size correction factor we applied to our 
measured values of r 50 % in order to recover actual values of r 5Q % 
for our objects. We apply a minimum and maximum correction 
factor of 1 and 2, respectively. 



3.3. Asymmetry 

Asymmetry is calculated by comparing an object with 
an image of itself rotated by 180°. It is therefore crucial 
to know the object's center, which becomes the pivot 
point. Our approach in determining centers m ostly fol- 
lows the one described in lConselice et all ((2000) . We first 
pick the brightest pixel, after some smoothing, as a first 
estimate of the center. We then refine it by calculating 
asymmetries within a 4-pixel radius aperture successively 
centered on each of the 9 points of a 3 x 3 grid at that ini- 
tial center. We proceed to the lowest asymmetry point, 
refining our mesh to sub-pixel level by interpolation until 
the difference in asymmetries between the lowest and sec- 
ond to lowest points is less than or equal to 0.001 or 20 it- 
eratio ns have been reached. Contrary to IConselice et al.l 
( 200(3) we only use a 4-pixel radius aperture, as opposed 
to full-aperture, to minimize our asymmetry, since our 
goal is really to find the center of symmetry of the bulge. 
Minimizing global asymmetry would generally give us a 
center closer to the center of light of the system, which, 
in the case of highly peculiar galaxies or mergers, could 
be very far from the bulge center. Although that ap- 
proach is just as valid, we found the first one to be a 
better discriminant of interacting systems. 

We implemented two other minor changes in the asym- 
metry algorithm. Both are mostly procedural, help re- 
duce the scatter, but also tend to produce values of 
A that are higher tha n that of the standard algorithm 
(|Conselice et al.| [2~000). though, we reckon, more accu- 
rate. One modification is the use of segmentation maps 
instead of circular apertures as mentioned above. More 
precisely, the procedure involves symmetrizing the maps 
first, and then applying them to the difference image ob- 
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tained by subtracting the rotated image from its original. 
This causes interacting systems to be fully included into 
the segmentation rather than only the part of which falls 
within a certain circular apperture. The second modi- 
fication we implemented conc erns the way we estim ate 
the effect of the background. (jConselice et al.l[2000t ) use 
a nearby empty region of space to calculate the asym- 
metry of the background and then subtract that asym- 
metry from the original value. We, on the other hand, 
estimate the effect of the background by calculating a 
second asymmetry value (A 1 ) from a convolved version 
of the object, with the following 5-point average convo- 
lution 



:(fi,j + fi+lj + fi-1,3 + fij + 1 + /ij'-l) (5) 



where fij represents the flux at the pixel of the 

image. If we assume that the intrinsic asymmetry of the 
light does not change in the weakly convolved version, 
and there is evidence for such an ass umption to hol d 
with even bigger convolution kernels (jConselicei l2003[ ) , 
the difference between the two asymmetry values must 
be entirely due to background, and since the standard 
deviation of the background in the smoothed image is 
reduced by a factor y/E from its unsmoothed version, we 
have 
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This implies that we can estimate the amount of asymme- 
try due to random fluctuations in the background from 
the following formula: 
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background 



A- A' 

1-1/75 



(9) 



and subtract it from our asymmetry measurement to ob- 
tain the intrinsic asymmetry of the object. This prevents 
us from subtracting background asymmetry from regions 
where it is due to intrinsic differences between opposite 
parts of the galaxy. 

3.4. dumpiness 

dumpiness is calculated by subtracting from an im- 
age a blurred version of itself. The blurred version is 
obtained by convolving the image with a circular tophat 
filter of diameter equal to 0.3 x r petro . After subtraction, 
only positive values are retained and summed. Straight 
forward application of this procedure almost always re- 
tains significant flux in the cente r, where bulges s harply 
peak, and must be corrected for. IConselicel (|2003h chose 
to simply blank the region inside one filter radius. We 
opted for a similar but slightly different approach. We 
decreased the size of our convolution kernel for points 
inside two filter radii. In that regime, we set the fil- 
ter radius to one half of the distance to the center. It 
reaches zero in the center, with the 9 central pixels not 
being smoothed at all, therefore always subtracting out. 
This allowed us to minimize the contribution of the cen- 
tral peak to the dumpiness value while still picking out 
bright clumps or bars near the center. In most cases 



though, this slight modification did not yield different 
results from the standard approach. As for the effect 
of background on du mpiness, we corr ected for it using 
the standard method (|Conselicel |2003[ ). We also make 
use, here, of segmentation maps for summing over galaxy 
pixels in the difference image. 

3.5. Gini 

The Gini coefficient ([Abraham et alj [2003) measures 
the inequality of the distribution of flux among the pix- 
els associated to a galaxy. Its possible range of values 
goes from 0, in the case where all the pixels would have 
the same intensity, i.e. complete equality among pixels, 
to 1, in the case where all the flux of a galaxy would 
be contained in a single pixel. In general, it can be de- 
fined as the ratio of the area between the Lorenz curve 
and the curve of uniform equality to the area under the 
curve of uniform equality. The Lorenz curve, L{p) 1 is in 
turn defined as the curve representing the proportion of 
the total flux contained in the dimmest p fraction of pix- 
els. The Gini coefficient is thus somewhat analogous to 
the concentration parameter, except that it is not mea- 
sured with respect to a specified center. It requires, on 
the other hand, objects to be segmented. In other words, 
a boundary needs to be drawn inside which pixels are as- 
signed to the obje ct. As described a bove (3.1), we follow 
the prescription of lLotz et al.l (|2004D and use the isophote 
of surface brightness fJ-(r petro ) in the convolved image as 
our boundary. We then calculated Gini coefficients for 
our objects by summing the values of the segmented pix- 
els in the following way: 
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where the n pixels are first sorted from dimmest to 
brightest (in absolute value), and Xi represents the flux 
of the i th pixel. This method is equivalent to the defini- 
tion given above, the absolute va lues making it fu rther 
more robust to background noise (jLotz et al.ll2004h . 

3.6. M20 

M20 stands for the normalized second order moment of 
the brightest 20% of the galaxy's flux. It is best described 
mathematically as: 



M 20 = log 10 



Si fi 



(11) 



where fi and represent the flux and distance from the 
center of the i th pixel respectively, and where the sum is 
performed by adding pixels in decreasing order of bright- 
ness (starting with the brightest one) until J2i fi reaches 
20% of the total flux. M tot in this equation is simply 
the second-order moment summed over all pixels. M20 
is thus like an inverse concentration for galaxies whose 
profile declines monotonically and isotropically. In those 
cases, the brightest 20% of the flux pixels is equivalent 
to the region enclosed by the 20% of the flux radius, and 
r 80 % and r petro then also follow a simple relation. How- 
ever, M20 is much more strongly influenced by bright 
clumps in the outskirts of galaxies than is concentra- 
tion. Therefore, whereas concentration can sometimes 
be thought as a bulge-to-disk ratio, M20 diverges from 
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Fig. 3. — The distribution of our full sample in concentration(C), asymmetry (A), clumpiness(S), Gini coefficient(G) and second moment 
of the brightest 20% of the light (M20). 



that concept in cases where extended non-axisymmetric 
light becomes important, as is the case in mergers and 
certain disks. 

4. RESULTS 

4.1. Morphology of z ~ 0.7 Galaxies 

Figure [3] shows the distribution of our full sample in 
our five morphological parameters. The distribution dis- 
plays a clear morphological bi modality, analogo us to the 
color bimodality observed by iBell et afl (|2004f ) at those 
rcdshifts. This bimodality appears in both asymmetry 
and dumpiness parameters, but becomes most apparent 
in a 2-D distribution when either of these is combined 
with a concentration-like parameter (concentration, Gini 
or M20). This implies that galaxies tend to be either 



spheroidal, in which case they have nearly zero asym- 
metry and dumpiness, or dominated by disks, which ex- 
hibit typical values of asymmetry and dumpiness around 
A = 0.25 and S = 0.17. Intermediate morphologies in 
which a visible disk still exists albeit being supplanted 
in importance by the bulge, such as in SO's and Sa's, are 
encountered less frequently. 

Although dumpiness and asymmetry appear to behave 
in very similar ways and even to actually correlate with 
one another, they are not completely degenerate. For ex- 
ample, the galaxies with the highest ratios of S/A are al- 
most exclusively edge-on galaxies, whereas the ones with 
the lowest ratios often have a bright compact center with 
a long but faint and smoot h tail, or cloud, exten ding on 
one side (tadpole galaxies [Griffiths et all Il994j would 
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Fig. 4. — Histograms of the relati ve number of ob jects in our sample as a function of eac h of our morphological parameters, compared 
to the distribution in the samples of Consclicc (2003) (light grey line) and Lotz et al. (2004) (dark grey line). 



fall in that category), thus boosting asymmetry, but not 
dumpiness. On the other hand, it is true that elliptical 
galaxies will have both near-zero asymmetry and dumpi- 
ness. 

The strong correlation between concentration and M20 
is due to the fact that they have similar definitions 
(see section 3.6). As mentioned in section 3.6, however, 
M20 is much more sensitive to non-axisymmetric features 
which occur in many disks as well as in mergers. This 
is why outliers lying above (higher M20) the M20-C re- 
lation begin to appear as one moves towards the low-C, 
high-M20 part of the plane. 

Figure [4] shows the distribution of our values of C, A, 
S, G and M20 comp ared to the one i n the samples of 
IConselicd (|2003[ ) and lLotz et "all (|2004l ). Because of the 
fact that we effectively use the flux within one petrosian 
radius in our estimate of the total flux ( see se ction 3.1) as 
opposed to 1.5 x r petro as in IConselicd ([2003] ) . our values 
of C are lower by about .3. From th eoretical consid- 
erations of Sersic profiles (|Sersid 1X9681 ). we obtained a 
similar shift rather ubiquitously across all Sersic indices 
going from 1 (exponential) to 4 (de Vaucouleurs). Other 
than that constant shift, all distributi ons are otherwise 
fairl y similar, except for the fact that IConselicel (|2003l ) 
and lLotz et al.l (|2004h have a higher fraction of objects 
with bulge-dominated morphologies. This, however, has 
no implications since their sample was selected by hand. 

The asymmetry values in our sample are also slightly 
different (AA ~ 0.1). This is certainly at least partly 
a real effect as high er redshift galaxies tend t o show 
more peculiarity (e.g. Ivan den Bergh et al.lll996h . How- 
ever, we did implement a different background asymme- 
try subtraction algorithm, which is specially written so 
that to avoid oversubtraction. Our clumpiness(S) mea- 
surements are also systematically lower. Such a decrease 
in S with redshift w as predicted and calculated, however, 
by IConselicd (|2003). and is a consequence of bright knots 
getting smeared out in images with lower resolution and 
lower signal-to-noise, the latter being a consequence of 
surface brightness dimming. In addition, the higher the 
real dumpiness, the larger this effect is. Our dumpiness 
values, thus rarely exceed 0.5. 

Our Gini values also tend to be lower, but we think 
this is also simply due to lower signal-to-noise. Figure^ 
shows the picture of NGC 4790, an Sd galaxy, taken from 



the Digitized Sky Survey. Because it has a very uniform 
disk, it's Gini coefficie nt is low, . 4 in t his case, which 
is the low-end limit in Lot z et ail (|2004j )'s sample, but 
which is not yet as low as many objects in our sam- 
ple. However, we found that by degrading its image to 
the level shown in figure [5Jd, we were able to reproduce 
an image with a Gini coefficient of 0.3, similar to the 
ones of figure [5k and [5ji which represent two G ~ 0.3 
galaxies from our sample. On the other hand, we find 
the Gini coefficient of ellipticals to be a fairly robust to 
resolution and signal-to-noise effects, which means that 
image degradation only stretches the low-end of the Gini- 
distribution, and indeed, this is what we observe in fig- 
ure |H Fin ally, our values o f M20 correspond very well 
to those of lLotz et al.1 ((2004) . The only difference being, 
as we mentioned earlier, their higher fraction of bulge- 
dominated objects. 

4.2. Physical Properties of z ~ 0.7 Galaxies 

Figures [6] and represent the distribution of sizes, 
masses, surface mass densities, star formation rates, spe- 
cific star formation rates and restframe FUV - g color 
in our UV-detected and full samples respectively. Size, 
mass and star formation rate measurements were per- 
formed as described in sections 3.1, 2.2 and 2.4 respec- 
tively. Surface mass densities were derived using the fol- 
lowing relation: 



where the factor of 0.5 accounts for the average effect 
of inclination. We also plot, in figure [5J the conditional 
plot of figure \7\ that is a plot where each column has 
been normalized separately, i.e. the value in each bin 
has been divided by the total number of objects in that 
same range of the independent variable. We also show 
in figure [5J for every column, the 10, 50 (or median) 
and 90% quantiles. These clearly trace out how the de- 
pendent variable varies specifically as a function of the 
independent variable as the effect of the number of ob- 
jects present at a given value of the independent variable 
is removed by the normalization. They also provide us 
with a sense for the spread in the relation. 

As our selection criteria are optimized for morpho- 
logical analysis, our sample is only complete for large 
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(c) (d) 

Fig. 5.— a) Sd galaxy NGC 4790. NGC 4790 has a Gini coefficient of 0.4 and is 19 Mpc away, b) Degraded image of NGC 4790. It's 
Gini coefficient is 0.3. c) Galaxy at z = 0.74 in the COSMOS field with a Gini coefficient of 0.3. d) Another galaxy from the COSMOS 
field with a Gini of 0.3. This one at z = 0.65. 



( r 50% si 3 kpc) and massive (logM, > 10) galaxies. 
We nevertheless observe the properties of galaxies in our 
sample to be, within our completeness limits, consistent 
with that of local galaxy samples in sizes, masses and 
surface mass densities. They do, however, have higher 
star formation rates. This increase in star formation 
rate also reflects itself in figure [9j which shows our star 
formation rate distribution in comparison t o the l ocal 
star formation rate function of iMartin et alj (|2005bf ). In 
fact, figure[9]demonstrates that our distribution is shifted 
by a factor of about 3.5 (0.55 in log space) which cor- 
responds to the overall increase in star formation rate 
density of the Universe between z — and z = 0.7 
( Schi minovich et"aTll2005l ). We demonstrate in the sec- 
ond half of this section that, by comparing to local sam- 
ples, this shift is ubiquitous among the full range of star- 
forming galaxies. As for the color distribution, it dis- 
plays bimodal ity (figure [3, simila r to that observed at 
low redshifts l|Strateva et al.l [200lh . as well as very dis- 
tinct red and blue s equences when plotted against mass 
(|Wyder et al.l l2006f ). with the transition occurring at 



logM* w 10.5 M . When plotted against star formation 
rate, it shows that the range of possible SFR's widens 
as one progresses to redder colors with the population 
bifurcating into old quiescent galaxies on one side, and 
dust-enshrouded star-forming ones on the other. Many 
of the relations in figures [6] t o [H have been studied in 
the low-redshift population b y lKauffmann et al.1 (|2003bl ) 
and lBrinchmann et al.l (|2004l ). We thus now turn to com- 
pare our results with theirs to investigate whether or how 
t hese relations have change d since redshift z ~ 0.7. 

iKauffmann et all (|2003bh observed log /i* to be propor- 
tional to log Af* for log M* < 10.5 followed by a flattening 
at higher masses. The relation between log //, and log M* 
in our UV-detected sample ( figure is consist e nt with 
theirs, but comparison of the [Kauffman n et all (|2003bl ) 
relation with our full sample (figures [7] and [5J is harder 
to reconcile. Within our detection limits, our two rela- 
tions are consistent at low-masses, but we fail to observe 
the break in the slope for higher-mass objects. However, 
because of our incompleteness at small radii, we are miss- 
ing the top part of the relation (where 750% ^1-0 kpc), 
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Fig. 6. — Distribution of our UV-detected sample in various physical parameters. Dotted lines represent in the logM, — log sSFR plane: 
our detection limit of log SFR = 0.11; in the log fi* — logr 5 Q% plane: our detection limit of log M, = 9.1 and the line of log M* = 11.6 
(which is the high-mass cutoff in our mass distribution); in the log/x* — logM* plane: our detection limit of r 50 % = 1.0 kpc; in the 
logM* — logr 5 Q% plane: the line of log/x* = 9.5 (which represents the typical value for ellipticals). Dashed lines represent in the log sSFR 
vs. logr 5 Q<% plane: the up per envelope of th e Brinch mann et all l)2004l 'l local relation shifted up by 0.45; in log sSFR vs. log M*: the upper 
and lo wer envelopes of the Brinchmann et al. (2004) local relation shifted up by 0.55; in log sSFR vs. log/i»: the upper and lower envelopes 
of the [Brinchmann et al. (2004) l ocal relation shifted up b y 0.45; in log SFR vs logM*: the median for blue galaxies as well as the upper 
envelope of the local relation of Brinchma nn et aLl ([2004) obtained from their raw fiber measurements (their figure 17), both sh i fted up 
by 0.35; and in log/x* vs logM*: the upper and lower envelopes as well as the median of the local relation from Kauffmann et al. (2003b) 
(unadjusted). 
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Fig. 7. — Same as figure [6] but for our full sample. Star formation rates below our detection limit have been obtained through SED 
extrapolation and are thus prone to large uncertainties. Dotted lines have been added in plots involving SFR that represent our detection 
limit of log SFR = 0.11. Addi t ional dashed-dotted lines are also shown that represent the (unadjusted) lower envelopes of the relations 
obtained by Brinchmann ct al. (2004) for local galaxies. 
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Fig. 8. — Conditional plot of figure [7] that is with each column being normalized separately. The solid lines represent the 10, 50, and 
90% quantiles in every column. We reiterate that star formation rates below our UV-detection limit have been obtained through SED 
extrapolation and are prone to large uncertainties. 



which deprives us of an upper envelope for the linear 
part of the relation thus making it hard to conclude one 
way or the other. A flattening at high masses would still 
be consistent with our results in the case that the upper 
envelope for galaxies at z ~ 0.7 be higher than in the 
iKauffmann et all (|2003bl ) relation. 

We also observe r 50 % to behave, as a function of mass, 
in a wa y similar to that repor ted by IKauffmann et al.l 
(2003b) for local galaxies and iBarden et al.l (|2005h for 
samples with redshift up to z ~ 1, namely slowly in- 
creasing at low masses followed by an accelerated in- 
crease and truncation of small size galaxies at masses 
> 10 10 ' 5 M with a typical spread of about an order of 
magnitude in r 50 % . The median of our relation is slightly 



shi fted towards large r radii c ompared to the trend found 
by IKauffmann et all (|2003bl ) for local populations, but 
this is simply an artifact of our incompleteness at low 
radii. Our d i stribu tion is close, though, to the one of 
IBarden et al.l (l2005[) for similar redshifts. 

Analogously to iBrinchmann et ahl (|2004f ) , we also find 
a correlation between log SFR and log M* at low masses, 
with a break at log M*/Mq > 10, where a bimodal i ty be- 
gins to emerge (figure |8|) . The IBrinchmann et ahl (|2004L 
figure 17) relation for fiber measurements of both SFR 
and M* augmented by 0.35 in log SFR is shown in dashed 
lines in figures [6l [7] and [8] (median and upper envelope 
representing 0.02 conditional probability contour) and is 
the relation that fits best our star-forming galaxies (UV 
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Fig. 9. — Star formation rate distribution in our UV-detected 
sample compared to the local star formation rate function of 
Martin ct al. (2005b) (solid line) and to the same function shift 
by 0.55 in log SFR (dashed line). 

sample). The dotted-dashed line shown in figures [7] and [8] 
represents the unmodified lower contour and appears to 
form a better envelope of our sample. It is thus possi- 
ble that only the star-forming galaxies have higher star 
formation rates at z ~ 0.7, though the star formation 
rates below our GALEX detection limit do carry large 
uncertainties. 

We do also observe the w e ll known, both lo- 
cally (|Brinchman n et ail l2004t |Percz-Gonzale z~et al.l 
20031) and at higher redshifts (iFeulner et all 120051 : 
Brinchmann fc Ellisl 120001: ICowie et all I1996D . relation 
between specific star formation rate and mass. Figures [7] 
and [8] again demonstrate a bimodality and spread in the 
relation at masses M*/M & > 10. We have seen that 
a shif t of 0.35 in log SFR from the Brinchmann et al. 
(|2004h fiber measurements fits best our data. They do 
calculate, however, an aperture correction (for both SFR 
and mass) and apply it to their study of specific star 
formation rate. Comparing with the aperture corrected 
data, th e shift necessary to bring t he log sSFR - log M re- 
lation of lBrinchmann et al.l (|2004 , figure 24) to fit ours is 
now 0.55 in log SFR, precisely the value representing the 
differenc e in psfr. of the Univers e between z ~ 0.7 and 
present (jSchiminovich et"aT1l2005h . Our specific star for- 
mation r ates are slight l y high er, however, than the ones 
found bv lFeulner et alj (|2005t ) for the same redshift. We, 
nevertheless, believe our results to be accurate as our 
distribution differs from theirs mostly at the high-mass 
end where we see a number of galaxies with significant 
(> 10 M Q yr -1 ) star formation rates, whereas their re- 
sults show an upper envelope at SFR 5 M Q yr^ 1 . As 
an aside, since specific star formation rate is correlated 
with the ratio of present to average past SFR, the be- 
havior we observe corroborates the evidence that more 
massive galaxies formed their stars, on average, earlier 
than lower-mass galaxies, also referred to as downsizing 



(jCowie et al.lll996t ). 

Our relation between specific star formation rate and 
mass surface density is also qualitatively similar to the 
local relation, with a transition at log/i» > 9 from 
star-forming to quiescent objects. This suggests that 
star formation ceases beyond a certain stellar den- 
sity, in agreement with [Kauffmaim et all (|2003b| ) and 
IBrinchmann et all (|2004D . More precise comparison fur- 
ther shows that our upper envelope for the distribution 
function is shifted up by 0.45 in log sSFR relative to 
IBrinchmann et all (|2004h 's as plotted in figures [6] to [H 
Similarly, the relation between specific star formation 
rate and size for our UV-detected sample is consistent 
with the local distribution upped by 0.45 in log sSFR 
(figure [6]). Although in both relations the lower contour 
appears not to change from the local distributions, the 
large uncertainties in the star formation rate estimates 
below our NUV detection limit make it hard to speculate 
on the evolution of low star-forming galaxies. 

Since our mass surface density and size distributions do 
not appear to evolve much from z = 0.7 to the present, it 
is safe to assume that the observed evolution in specific 
star formation rate simply reflects the evolution in star 
formation rate itself. It thus appears that the relations 
between physical parameters of z ~ 0.7 galaxies change 
little relative to local samples but for an overall increase 
by a factor of 10 35 to 10 55 in their star formation rate. 
As lO 55 precisely represents the decline in star forma- 
tion rate density in the Universe from z — 0.7 to z — 0, 
this overall dimming with time appears to account for 
most of the star formation rate dens ity evolution sinc e 
that redshift, confirming the results of lWolf et alj (|2005l) . 
On the other hand, since our observed shift could be as 
low as 10 35 , our results can easily accomodate a mild 
number evolution in the form of a steepening of the faint- 
end slope in the FUV-l u mino sity function such as ob- 
served bv lArnouts et al.l (|2005f ). 

4.3. Relation of Morphological Parameters to Physical 
Properties 

In this section we look at how morphological param- 
eters relate to physical properties of galaxies. The 
most important result we discover is the fact that all 
morphological parameters show a transition at masses 
of 10 10 5 M Q and above, which we interpret as a shift 
from disk-dominated to bulge-dominated galaxies. This 
transition mass i s ident ical to the one observed by 
iKauffmann et al.1 (|2003bh for local galax i es and consis- 
tent with that observed by iBundv et alj (|2005h also in 
the redshift range 0.55 < z < 0.8. We also discuss in this 
section how incompleteness in the physical properties dis- 
cussed in the previous section affects our distribution of 
morphological parameters. 

Figures [10] and [TT] display each of our morphologi- 
cal parameters as a function of size, mass and surface 
mass density for our full sample, while figure [T^l com- 
pares the distribution in morphology of galaxies from our 
UV-detected and non-UV-detected samples. The large 
scatter in all the relations shown in figure \W\ indicates 
that no one morphological parameter is a good indica- 
tor of the physical state of a galaxy. Nevertheless, we 
do observe some broad correlations as well as a clear bi- 
modality in some cases, both of which we discuss in this 
section. We focus here on figures ITOl to [T2l and postpone 
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Fig. 10. — Morphological parameters as a function of radius, 
mass and surface mass density for our full sample. 



discussion of how star formation relates to morphology 
to the next section. 

Concentration appears to show some kind of relation 
to most of our studied physical parameters, the most in- 
teresting one being that with stellar mass. Figures [10] 
and [Tl] indeed show how the distribution of galaxies 
in the C — log M* plane displays a sharp transition at 
log M*/M© = 10. 5. This is much t he sam e behavior 
as seen locally by iKauffmann et al.l ()2003bl ). and indi- 
cates a shift from a disk-dominated population for which 
C ~ 2.6 independently of mass to a bulge-dominated 
population whose concentration rapidly increases with 
mass. We find C to also behave in a similar way with 
log/z*. There appears, on the other hand, to be little 
correlation between concentration and size in our sam- 
ple for the most part except for the fact that the point- 
spread function has a greater effect on the central profile 
of smaller galaxies, which causes the observed cutoff at 
progressively lower concentrations as we go to smaller 

^50%- 

We have already discussed how the distributions of 
asymmetry and dumpiness show a bimodality with 
bulges and disks segregating at the low and high-end re- 
spectively. This can be seen again in figure [12] as well 
as in most relations of A and S in figure [TO] For exam- 
ple, although asymmetry and dumpiness do not seem to 
correlate with stellar mass in neither disk-dominated nor 
bulge-dominated galaxies, the two populations differen- 
tiate in their respective values of A and S, bulges lying 




0.0 0.2 0.4 0.6 0.8 10 11 8 9 
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Fig. 11. — Conditional plot of figure flOl with each column being 
normalized separately and lines representing the 10, 50, and 90% 
quantiles in every column. 



almost exlcusively in the range of to 0.1 for A and -0.05 
to 0.05 for S, and disks spanning almost the entire range 
of possible values, with typical values around A = 0.25 
and S = 0.17. We again see the transition from a disk- 
dominated population to a bulge-dominated one to take 
place at masses of around 10 10 ' 5 M Q . Although a bi- 
modality can also be seen in the A or S vs. log /i» planes 
(figure [TO]), the two populations definitely follow a more 
linear relation with both A and S monotonically declin- 
ing with log /i» (figure Hip . Asymmetry and dumpiness 
also appear to correlate with size in disk galaxies. We 
demonstrate in the appendix how such a relation between 
dumpiness and size is expected in a resolution-limited 
sample. As for asymmetry, it is unclear whether the ob- 
served trend is real or artificial. Certainly, our brightness 
cut of / < 23 mag corresponds, at any given radius, to 
a surface brightness limit which progressively increases 
as we go to smaller radii. As spheroids have higher sur- 
face brightnesses than disks, this selection effect can skew 
the observed population at a given size, which in turn in- 
fluences the median asymmetry, and contributes to the 
observed trend. 

Beyond the bimodality, we also note the crest extend- 
ing to high values of A in all the plots of A versus phys- 
ical parameters. This crest is interesting since it repre- 
sents galaxies that depart from the standard disk model 
by displaying asymmetric features, whether intrinsic or 
due to interactions. It lies straight above the normal 
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Fig. 12. — Distribution of our sample in our five morphological parameters. Black histograms show the distribution of the entire sample, 
while dark and light grey histograms respectively represent objects detected and objects not detected in the UV. 



disk population, implying that these galaxies do not have 
physical properties different from those of regular disks 
(though other properties, that we do not study, such as 
IR-luminosity or dynamics, would probably differ). 

The Gini coefficient is the only morphological param- 
eter to correlate with mass throughout our entire mass 
range, even though the steepest increase still occurs at 
masses between 10 10 5 and 1O U M0. On the other hand, 
Gini is mostly insensitive to the size of the galaxy, and 
hence correlates equally well with surface mass density. 
Nonetheless, we observe a lack of low- Gini small objects 
which is due to surface brightness incompleteness com- 
bined with the fact that the Gini coefficient systemati- 
cally increases at low resolution where objects approach a 
point-like psf profile. We also observe an increase, with 
size, of the low-end cutoff of the Gini distribution at 
larger radii (logr 50 % > 0.4). This increase appears to 
be real, and signifies that larger disks are allowed less 
uniformity in their light distribution, implying they are 
either more bulgy, more clumpy or have, in general, more 
structures, such as spiral arms. Often, all three apply. 

Moving to M20, as shown in figure [3j it is largely anti- 
correlated with concentration. It hence displays many 
of the same trends as concentration, only mirrored. For 
example, M20, just like concentration, appears to be in- 
dependent of mass for log M* < 10.5, but strongly (anti- 
)correlated with it at higher masses. A notable differ- 
ence, however, is that, as seen in figure [T2l star-forming 
galaxies occupy two thirds of the possible range of values 
of M20 (-2.0 to -0.5 relative to a full range going from 
-2.75 to -0.5), whereas they occupy only half of the pos- 
sible range of values of concentration (from 2.0 to 3.25 
relative to a range of 2.0 to 4.5). This is because M20 
is more sensitive than C to bright features far from the 
center. This reflects in the relation between M20 and 
logr 50 %, which shows a bifurcation. As discussed pre- 
viously for concentration, the downward trend in bulge- 
dominated galaxies is simply due to our resolution limit. 
On the other hand, just like in the case of asymmetry, it 
is unclear whether the upward trend in disk-dominated 
galaxies is real or also a consequence of our selection cri- 
teria. For one, the smallest diffuse disks do not make 
our magnitude cut, which eliminates high-M20 objects 
at low-r5o% , and then M20 and dumpiness are somewhat 
correlated so that resolution effects could also help cre- 
ate the observed trend. Nevertheless, though it might be 
accentuated by selection criteria at small radii, it is likely 
that the relation be real as it holds at logr 50 % > 0.5, our 



completeness limit, which suggests that larger disks are 
more likely to have bright regions in their outskirts, and 
that smaller galaxies have the bulk of their light distri- 
bution more centralized. 

In summary, we have seen in this section how all mor- 
phological parameters show, as a function of stellar mass, 
a transition at M* ~ 1O 1O 5 M0 which we interpret as a 
transition from disk-dominated to bulge-dominated ob- 
jects. We have also shown how the Gini coefficient is the 
only morphological parameter to correlate with stellar 
mass throughout our mass range. 

4.4. Morphology of Star-forming Galaxies 

4.4.1. Star Formation Rate as Related to Morphological 
Parameters 

The relation between color, star formation rate and 
specific star formation rate, and morphological parame- 
ters is shown in figures fT3l and [T4l for galaxies detected in 
the UV. As expected, disk galaxies are blue and become, 
on average, redder as the bulge becomes more prominent, 
that is as we move towards bulge-like morphologies. Nev- 
ertheless, blue galaxies seem to exist for all values of our 
morphological parameters, including a population of blue 
compact galaxies. Since many morphological parameters 
suffer from incompleteness at low radii, we also plotted 
in figures [15] and Ql)] the same relations but for galax- 
ies with logr 50 % > 0.55 (or r 50 % > 3.55 kpc) only, and 
still taken from the UV-detected sample. This cut rep- 
resents our size completeness limit. Relations of color to 
morphological parameters are stronger in this complete 
sample, but their qualitative behavior remains the same. 
One notable difference is the absence of blue compact 
galaxies, due to the simple fact that they do not make 
the size cut. 

From figures [T4l and [TBI it is clear that, among our five 
morphological parameters, the strongest relation color 
possesses is that with the Gini coefficient. Color is also 
strongly correlated to asymmetry as well as to dumpi- 
ness for A < 0.2 and S < 0.12. However, it then quickly 
stabilizes and changes very little for values of A > 0.4 and 
S > 0.2. These strong correlations at lo w values of A or 
S resp ectiv ely confirm simi lar results of IConselice et all 
(2000) and IConselice] ([2003'). The upper ranges of asym- 
metry and dumpiness, however, have not been probed 
before. We find, here, that those highly asymmetric or 
very clumpy galaxies possess the same blue colors as their 
more standard disk counterparts. Similarly, we find no 
significant difference, as discussed below, in the star for- 
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Fig. 13. — Absolute and specific star formation rates as well as FUV - g color of our UV-detected sample as a function of morphological 
parameters. 



mation rate or in the specific star formation rate of high- 
A or high-S galaxies compared to normal disks. We dis- 
cuss our interpretation of the relation between color and 
morphology further in section 5. 

Interestingly, the parameter with which star formation 
rate correlates the best is the Gini coefficient. Figure H4l 
indeed shows that the logarithm of the star formation 
rate is a linear function of G (with log SFR oc 2.3 x G), 
although the relation has a dispersion of as much as 1 
to 2 orders of magnitude (at the 80% confidence level) 
from the lower to the upper end of the range in G. In 
general, quiescent early-type galaxies have higher values 
of G than most of their star-forming counterparts, and 
actually dominate the high-G population as shown in 
figure [T2l Star- forming galaxies that do, however, attain 
equally high values of G, seem to typically carry star 
formation rates > 10 M© yr _1 , which suggests that the 
strongest starbursts have very luminous components that 
stand out from the rest of the galaxy. From inspection, 
we notice that these are often central components, pos- 
sibly in the process of being formed. We thus interpret 
this result as a signature that episodes of strong star for- 
mation can be often linked to bulge growth, growth that 
sometimes turns galaxies into early-type objects. We, 
again, address this subject in more detail in section 5. 

Otherwise, the star formation rates for our UV- 
detected galaxies, do not show a significant correlation 
with any other one of our morphological parameters, in- 
cluding asymmetry and dumpiness (though a slow de- 
cline of SFR with M20 is visible in our size-complete 
sample). Given that all of our morphological parame- 



ters are normalized quantities, one might expect specific 
star formation rate, which is also a normalized quan- 
tity, to show a stronger correlation with morphologi- 
cal parameters, and in particular with dumpiness which 
wa s shown t o corre late with the Ha equivalent width 
by iConselicel ((2003) . This is, however, hardly the case. 
Specific star formation rate does show lower values for 
bulge-like morphologies and a rapid rise towards disk- 
like morphologies, but it then quickly tappers off so that 
among disks, the sSFR becomes rather independent of 
mo rphology. For en ds of comparison, we have also plot- 
ted [Conselicc (2003)'s relation between sSFR (converted 
from FiW(Ha) and renormalized to redshift z ~ 0.7 by a 
factor of lO 35 ) and dumpiness, in figures [15] and [TCI 
This shows that our two results are not inconsistent, 
only that there exists a larger variety of gala xies, espe- 
cially at z <~ 0.7, than was first considered by IConselicel 
(2003). In particular, these include large spirals which, 
although harboring substantial star formation, tend to 
have lower specific star formation rates because of their 
large masses (cf. the sSFR-mass relation in figure [7]) . 
Their spiral structure and star formation regions would, 
nevertheless, confer them high values of dumpiness. On 
the contrary, lighter, more compact and more floccu- 
lent late- type disks, would typically have higher specific 
star formation rates than the typical grand-design spiral, 
though often lower values of dumpiness. This diversity 
in the galaxy population adds scatter, so that it becomes 
difficult to infer specific star formation rates from mor- 
phology alone. We further discuss the relation between 
specific star formation rate and dumpiness in the ap- 
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Fig. 14. — Conditional plot (i.e. with each column separately normalized) of the absolute star formation rate, specific star formation 
rate and restframe FUV - g color of our UV-detected sample as a function of morphological parameters. The lines represent the 10, 50 and 
90% quantiles in every column. 



pendix. 

4.4.2. Morphology, Color and Specific Star Formation Rate 

In the last two sections we have investigated relations 
between morphological and physical parameters, and al- 
though we discovered some interesting trends, inference 
of physical properties of galaxies from any one morpho- 
logical parameter is uncertain at best, because of the im- 
portant scatter in all of the existing correlations. We now 
examine whether the situation can be improved upon 
by considering two morphological parameters simultane- 
ously We have already demonstrated in section 4.1 how 
this can help separate spirals and ellipticals. We now 
take one step further, and study color, specific star for- 
mation rate and FUV attenuation as a function of multi- 
parametric morphology. 

Our approach has been to color code our morphology- 
morphology graph (figure[3|) with the median of the stud- 
ied quantity (FUV - g color, sSFR or Apuv) in a number 
of bins across the distribution. The results are shown in 
figures [17] through [19] On the diagonal, we show the 
median distribution of FUV - g color, sSFR and Apuv 
successively as a function every one morphological pa- 
rameter. We note that bins with fewer than 10 objects 
were discarded. 

Figure [T7] demonstrates how, in the main galaxy pop- 
ulation, the color and morphology bimodalities strongly 
overlap. By comparing figure \T7\ to figure [3] one can see 
indeed, that the regions of morphological space where 
disk-dominated galaxies lie are blue and the ones pop- 



ulated by bulge-dominated ones are red. Although this 
relation between mor phology and color has been known 
for a long time (e.g. Ide Vaucouleursl fl961h . we hereby 
demonstrate the strength of that correspondence in the 
particular case of FUV - g, as well as in the sharpness of 
the transition between the two populations. Indeed, in 
most graphs of figure 1171 the median color jumps from 
FUV - g ~ 2 (blue) to FUV - g ~ 5 (red) over only 2 
or 3 morphological bins. It is otherwise very little cor- 
related with morphology within the two distinct popula- 
tions themselves, and the color is thus very uniform over 
both regions of disk (blue) and bulge-dominated galaxies 
(red). We note, however, that almost all of our bins, in- 
cluding the bulgy morphologies, do contain at least some 
blue galaxies. In the bulge-dominated bins, these are the 
compact blue galaxies mentioned earlier. Although very 
interesting, they are far less numerous and the majority 
of galaxies follow the pattern described above. As a fi- 
nal point, we mention, that the converse, which would 
be bins of disk-dominated morphologies containing red 
galaxies, is not seen (figure |2H|) . 

Although most graphs in figure [T7] show quite a sharp 
transition, color seems to correlate more smoothly with 
the Gini coefficient than with any other morphological 
parameters, so that graphs that include G, but especially 
the A-G and S-G planes, show wider transition regions. 
In the S-C and S-M20 planes, the region of concentrated 
galaxies that possess high dumpiness also shows inter- 
mediate median colors. By inspecting the distribution of 
colors in those transition bins (as an example, we show 
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Fig. 15. — Absolute star formation rate, specific star formation rate and restframe FUV - g color of our UV-detected size-complete 
sample as a function of morphological parameters. The dashed line in sSFR vs. S represent the relation of Consclicc (2003) normalized to 
2 = 0.7 and to our distribution in S. 



the color distribution in bins of the A-G plane in fig- 
ure [20]), we discover that it is very wide, with all colors 
from FUV - g = 1 to FUV - g = 6 being represented more 
or less equally. This indeed results into a typical median 
around 3.5, i.e. green. As we move to bluer colors, we see 
that we begin to progressively lose the reddest objects, 
until we are left with only blue galaxies (—0.5 < FUV 
- g < 3) in the dark blue bins. In the red bins, on the 
other hand, red galaxies dominate over the other ones, 
though as just discussed, some blue galaxies might still 
be present. 

The graph of specific star formation rate as a function 
of morphology (figure [T8f is qualitatively similar to the 
color-morphology one (figure [T7|) . However, its interpre- 
tation is a little different since the sSFR distribution is 
substantially more bimodal than the color distribution 
(figure [7]). In that context, the median is more indica- 
tive of the number ratio between the two populations 
rather than the typical specific star formation rate of a 
galaxy in that bin. The mode of the sSFR distribution 
of the quiescent population is situated, in our poor es- 
timate, at around log sSFR = —11.5 (figure [7]), value 
that is represented by red in figure [18] This means that 
bins that include only quiescent galaxies are going to be 
red. Orange bins, on the other hand, are likely already 
to contain a few star-forming objects, yellow bins even 
more, and so on. The mid-point, where there is about an 
equal share of star-forming and quiescent galaxies, lies at 
log sSFR ~ —10.4: green. Once one reaches dark blue 
bins (log sSFR > —9.3), only star-forming galaxies re- 



main. By looking at the various plots in figure 1181 it 
thus becomes apparent that most regions of morphology 
space, including the bulge-dominated areas, do contain 
at least some fraction of star-forming galaxies. Some 
of the star-forming galaxies with bulge-like morphologies 
are the blue compact objects, but many of them are, in- 
stead, red and dusty. This is demonstrated in figure HH 
where we show the FUV— g color of galaxies possessing 
a Gini coefficient of G > 0.57, as a function of their star 
formation rate. It also shows that the red and dusty 
galaxies are the ones that have the highest intrinsic star 
formation rates in our sample, though because they also 
tend to be massive, their sSFR is still typically lower 
than that of less massive galaxies (figure [7J . 

To confirm the presence of dusty star-forming galaxies 
in bulge-dominated regions of morphology space, we plot- 
ted in figure [H] the relation of Apjjy to morphological 
parameters for our UV-detected sample. All graphs in- 
deed show a rise of Afuv in those regions. These trends 
seem to indicate that the amount of attenuation could 
be mostly due to geometrical effects: compact galaxies 
(high C and G, and low A, S and M20 galaxies) hav- 
ing the highest extinction and galaxies with extended 
star-forming disks or interacting objects with disturbed 
morphologies (high M20 or high A objects) possessing 
lower values of Apuv ■ However, as the regions of high 
extinction seem to coincide with those where spheroids 
are located, it is possible that they be contaminated by a 
number of bulge-dominated objects mimicking the effect 
of dust with older populations and thus modifying the 
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Fig. 16. — Conditional plot (i.e. with each column separately normalized) of color, absolute and specific star formation rates of our 
UV-detected and size-complete sample as a function of morphological paramet ers. The s olid lines represent the 10, 50 and 90% quantiles 
in every column. The dashed line in sSFR vs. S represent the relation of Consclicc (2003) normalized to z = 0.7 and to our distribution in 
S. 



observed relation between attenuation and morphology. 
From inspection, we find that 77 out of 189 or ~ 40% of 
the UV-detected galaxies with FUV - NUV > 1 (which 
translates to Apuv > 4.12) in our sample have early- type 
morphologies; the remaining 60% show dusty-looking 
disks, though often also accompanied by an important 
central component. Despite the fact that this is poten- 
tially substantial contamination, we do not believe this 
effect to significantly alter the trends in figure [19] since 
dusty objects often possess larger UV-slopes than old 
populations anyway. In order to confirm the robustness 
of our observations, we tried to identify and remove ellip- 
tical galaxies present in our UV-detected sample by, first, 
selecting them through various cuts in the space spanned 
by our morphological parameters, and then through vi- 
sual inspection, but we always arrived at results very 
similar to those of figure 1191 

Dusty starburs ts are often luminous in infrared light. 
lLotz et ail (|2004T ) have shown that the most extreme of 
these IR-luminous galaxies, ULIRGs, typically have high 
R-band values of M20, as they often carry signatures of 
interactions, but also high values of G, which is indicative 
of the fact that their optical light primarily comes from 
one or a few regions of concentrated luminosity (such as a 
nuclear starburst). In our sample, the upper-right ridge 
of our M20-G plane shows higher FUV attenuation (fig- 
ure [in]) suggesting that galaxies in that part of the plane 
tend to be dustier. Because of dust, galaxies that carry 
such morphologies also display redder colors than typical 



blue-sequence galaxies (figure H~7|) even though they pos- 
sess similar specific star formation rates (figure [18]) . Our 
results thus appear to be consistent. We also observe 
that many strongly star-forming galaxies lie in the same 
regions of morphology space as bulges. These most likely 
represent more advanced stages of merging in which the 
resulting galaxy is finalizing its formation process. Alter- 
natively, these could also have grown from secular insta- 
bilities. We explore the nature of starbursting galaxies 
in more detail in section 5.2. 

It thus appears that the regions describing spheroidal 
galaxies (high-C, high-G and low-A, low-S and low- 
M20) are actually populated by a variety of objects: 
dusty red starbursts, blue compact galaxies, quies- 
cent elliptical galaxies, as well as UV-detected galax- 
ies that, from visual inspection, also appear to be el- 
lipticals. We consider three possible origins for the 
presence of UV-lig ht in elliptical galaxies: residual 
star formation (e g. iTeplitz et all [20061: lYi et all 120051: 
Salim et all 1200 5: Stan ford et all 2004; M enanteau et all 



20011 ). FUV-uptu r n from evolved hot horizontal-branch 



stars (lO'Connelll Il999t iBrown et alj 120001 : iRich et all 
[20051: iBoselli et all 120051 ), or weak AGNs. We immedi- 
ately rule out the FUV-upturn since it had been shown to 
fade rapidly with redshift reaching colors of FUV ' —V = 7 
at redshifts of z > 0.3 (|Lee et al.l 120051 and references 
therein), implying that we could only detect the FUV- 
upturn in galaxies with mi < 18.5 which is brighter 
than any galaxy in our sample. We cannot rule out 
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Fig. 17. — Restframe FUV - g color as a function of each as well as of each combination of two morphological parameters for 
our full sample. The value in each color bin as well as in each histogram bin was obtained by taking the median FUV - g color 
of all the objects present in that bin. Bins with less than 10 objects were removed. The color scale for histograms is indicated 
by the thin font top-most y-axis on the right side of the figure as well as the bottom-most y-axis on the left side of the figure. 
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Same as figure \T7\ but weighted by specific star formation rate. 
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Fig. 19. — FUV attenuation as a function of morphological parameters for UV-detected objects. Redder represents more 
attenuation while bluer represents less dust correction. The color displayed corresponds to the median value in that bin. Bins 
with less than 10 objects were removed. Histograms similarly represent median Afuv values across bins of each of our single 
morphological parameters, and are plotted on the scale indicated by the thin font top-most y-axis on the right side of the figure 
and the bottom-most y-axis on the left side of the figure. 
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Fig. 20. — Rcstframc FUV - g color distribution of our full sample in bins of the A-G plane. Every histogram spans a color range of 
-1 < FUV - g < 8. 
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Fig. 21.— FUV-g color versus log SFR for UV-detected , 
with a Gini coefficient of G > 0.57. 
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any of the remaining two options, but we do see signs, 
in the asymmetry and/or dumpiness residual images, of 
faint underlying disks in about 50% of our UV-detected 
spheroids with FUV — NUV > 1, suggesting that resid- 
ual star formation does play an important role. Al- 
though these galaxies appear to populate the lower half 
of star formation rates among UV-detected high- Gini ob- 
jects (the upper half being mostly populated by dusty 
nuclear starbursts), their dust attenuation (and hence 
their SFR) could still be overestimated, particularly for 
galaxies with significant 500 to 800 Myr old popula- 
tions, since such populati ons can carry steep UV-slopcs 
(|Bruzual fc Charlo 



opuiatK 
tj l2003ft 



4.4.3. Stellar Mass-limited Sample 

As many studies nowadays utilize stellar mass-limited 
samples, we investigated how a mass cut would affect our 
results, focusing on figure [HI the graph of specific star 
formation rate as a function of morphology. We chose 
a stellar mass cut of log > 10, which is roughly our 
mass completeness limit. Only in two instances do we 
observe significant differences with figure [T5J in C vs G 
and in M20 vs G. The two plots for this mass-limited 
sample are shown in figure [22] along with the histogram 
of the median sSFR as a function of G. We also show 
in figure [221 the same plots again for our full sample in 
order to facilitate comparison. We observe that, in our 
mass-limited sample, the median sSFR flattens for val- 
ues of G < 0.5, whereas it keeps rising at low values 
of G in our full sample. The shape of the contours in- 
dicate that this difference originates primarily from the 
regions of the plots where C » 3 or M20 » —2 and 
0.4 < G < 0.55. We think that such a difference in 
behavior at intermediate to low values of G can be sim- 
ply explained by the fact that lower mass objects have 
lower signal-to-noise, and that galaxies with low signal- 
to-noise tend to have lower Gini coefficients (cf. §4.1). 



Lower mass galaxies thus tend to populate the upper-left 
half of the C-G plane, or the lower-left part of the M20- 
G plane, in higher numbers compared to the lower-right 
and upper-right halves respectively. Because they have, 
on average, higher specific star formation rates than more 
massive galaxies (figure [7|), they tend to augment median 
star formation rates in those regions favorably relative to 
other parts of the planes, resulting in the behavior seen 
in figure 1181 In any event, as these differences between 
our full and mass-limited samples remain small, all the 
important results we have discussed so far apply equally 
to both samples. 

5. DISCUSSION 

5.1. Implications for Blue to Red Sequence Evolution 

The results discussed in the previous section point us 
towards a certain scenario of galaxy evolution. The 
strong correspondence between color and morphology 
suggest that the two characteristics evolve concurrently 
as we, indeed, observe no red disk-dominated galaxies 
and few bulge-dominated blue ones. The scenario is more 
ambiguous, however, when one looks at the star forma- 
tion rate, absolute or specific, as a function of morphol- 
ogy. This is due to the fact that, as we have shown, there 
exists a segregation in the star formation rate of bulge- 
like galaxies with, on one hand, dusty nuclear starbursts 
(blue compact galaxies can also have high star formation 
rates), and, on the other hand, old quiescent bulges. This 
split is not present in color since, except for the blue com- 
pact galaxies, all of the above types of objects are red, 
which is why color shows a better correspondence with 
morphology than does specific star formation rate. This 
implies that neither color nor morphology alone are good 
indicators of the physical state of a system. Instead, ge- 
ometry and dust conspire to link color and morphology 
in a way that appears independent of the star formation 
history. 

The fact that many among the strongest star-forming 
galaxies possess bulge-like or near bulge-like morpholo- 
gies, as measured by their Gini coefficient, indicates that 
episodes of strong star burs t s are often linked to bulge 
growth. iMihos fc Hernquistl (|1996f ) showed that gas-rich 
major mergers, that can tra nsform disk galaxies int o el- 
lipticals, trigger star bursts. iRobertson et al.l (|2006j ) fur- 
ther showed that even when the gas is not entirely con- 
sumed and that a disk remains at the end of the pro- 
cess, considerable bulge growth still occurs. It is thus 
likely that many of our strongly star-forming galaxies 
are late-stage mergers, and red-sequence objects in the 
becoming. The fact that we also observe many of these 
high-SFR objects to be red i s consistent with them of - 
ten being IR-lum i nous ( e.g. ISanders fc Mirabel 1 1996t ). 
IChakrabarti et all (|2006f) indeed showed that gas-rich 
mergers emit most of their bolomctric luminosity in the 
IR after the first pass due to infall of gas towards the 
center causing increased extinction. Such merger events 
therefore appear to preserve the color-morphology corre- 
spondence mentioned above, and could explain our ob- 
servations. 

Semi-analytical models (|De Lucia et al.l 120061 ) indi- 
cate that some low to intermediate-mass field ellipticals 
(M* < 10 n M Q ) could have built a non-negligible frac- 
tion of their mass from a last major merger event at 
redshifts as low as z = 0.5. This is thus consistent with 
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Fig. 22. — Median specific star formation rate as a function of C and G (first column), M20 and G (second column), and G only (third 
column), for galaxies with log Af* > 1 (to p row) as compared to the same relations in our full sample (bottom row). Note that the bottom 
row plots are identical to their figure [T81 counterparts. We show them again here to facilitate comparison with our mass-limited sample. 
Also note the different axis scales on the two histograms. 



such objects b eing present in our sample. On the obser- 
vational front. [Bund v et all (|2005T ) found that, although 
the integrated stellar mass density evolved only mildly 
since z ~ 1, its morphological mix underwent significant 
changes, also hinting towards a significant occurrence of 
merging since that redshift. Studies of B an d R-band lu- 
mino sity functions deri ved from the DEEP2 (iFaber et al.l 
120051) and COMBO-17 (|Bell et al.ll20"ollWolf et al.ll2003h 
surveys similarly point towards such a scenario. 

Others, on the other hand, find evidence that the most 
massive galaxies were already in place by redshift z ~ 1 
(I Juneau et all 2005: Tre u et al.|[2005HFeulner et al.ll2005t 
iBrinchmann fc Ellis 20Q3). We, on the contrary, observe 
that ~ 30% of our objects with mass log M„/Mq > 
10.5 are experiencing significant star formation (SFR 
> 1 M Q yr _1 ), in close agreement with lBell et al.l (|2005h . 
and in support of the scenario in which a fraction of red 
sequence galaxies continue to grow after z ~ 1. 



5.2. Nature of Starbursting and Blue Compact Galaxies 



In order to investigate the nature of the starbursting 
galaxies in our sample, we pulled out all galaxies with 
a SFR > 30 M yr _1 , and performed a visual classifi- 
cation on them. We find ~ 30% of them to be normal 
spirals, 25% appear to be ongoing a major merger event, 
while another ~ 35% have bulge-like morphologies with 
either a clear underlying disk or features indicating re- 
cent events. The remaining 10% of our high SFR ob- 
jects appear to possess either compact or regular ellip- 
tical morphologies, and could be AGN contaminations, 
or simply nuclear starbursts without any morphological 
peculiarities. 

As mentioned in the previous section, strongly star- 
forming galaxies with bulge-like morphologies could be 
objects in late stages of a merger where only one newly 
formed galaxy remains with a strong nuclear starburst, 
a large bulge and most debris having either fallen back 
onto the galaxy or having been dispersed. Alterna- 
tively, it is also possible that some of these bulges are 
being forme d by other means suc h as secular insta- 
bilities (e.g. iDebattista et all 120061 : IWang et all 120061 : 
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iKormendv &: Kennicutdl2004h . also suggested by the fact 
that most of our star bursting late-type spirals possess 
bars. Nonetheless, the existence of these objects ap- 
pears robust as they h ave also be e n obse rved in infrared 
studies. For example, iBell et ail (|2005l ) classify a frac- 
tion o f their strongest starbursts into the E/SO cate- 
gory. [Melbourne et al.l (|2005( ) similarly finds a signifi- 
cant fraction (~ 20%) of luminous infrared galaxies at 
0.61 < z < 1. 00 to posses s comp act or elliptical mor- 
phologies, and IZheng et all (|2004l ) finds ~ 25% of their 
LIRGs to have compact morphologies. Although it is 
difficult to compare our UV-based measurements with 
these numbers, they are still mostly consistent. The 
presence of such objects has implications on star forma- 
tion quenching time-scales, and might p ose some contra- 
dictio n with the merger simulations of iDi Matteo et aTl 
( 2005) in which the gas expulsion and star formation cut- 
off occurs before the two galaxies fully come together and 
settle into a spheroidal system. 

We also mentioned the existence in our sample of a 
population of blue compact galaxies, a term we have 
loosely used to describe blue galaxies (FUV - g < 2) 
with either high concentration or Gini, or low asym- 
metry, dumpiness or M20. Although this is a some- 
what different defin ition from the one usually used in 
the literature (e.g. iNoeske et al.l 120061 and references 
therein) which is based, among other things, on size 
(typically r 50 % < 3.5 kpc) and surface brightness (typi- 
cally ft(< r 50 %) < 21 B mag arcsec -2 ), most objects we 
characterize as blue compact galaxies do possess half- 
light radii and surface brightnesses that fit the above 
criteria. These objects are thought to be dwarf galaxies 
brigh t ened by intense star f ormation (e.g. INoeske et al] 
[20061: IGuzman et all 12001 IGuzman et all 119971) . A 
few blue concentrated galaxies, however, do show ex- 
te nded componen t s, and cou ld represent, as suggeste d 
bv lHammer et al.l (|200lD and lBarton fc van Zed ([2001h . 
bulges of spirals, in formation. We cannot rule out, 
however, the possibility that some of them harbor weak 
AGNs which would confer them their blue FUV - g color. 
This question will be further investigated in future work. 

5.3. How Independent are Measurements of Star 
Formation Rate and Morphology? 

Throughout this paper we have attempted to look at 
how star formation rate and galaxy morphology corre- 
late. The question can be asked however as to whether 
these are truly independent quantities, or simply two 
measurements of the same phenomenon. Is the light dis- 
tribution representative of the structure of the galaxy? 
Evidently, part of the answer is no, as star formation 
regions tend to stand out on top of the underlying struc- 
ture. This effect becomes all the more prominent as one 
moves to bluer bands. In addition, significant structural 
components of a galaxy might be hidden under low sur- 
face brightness, and elude detection. 

Clumpiness, for one thing, is precisely an attempt to 
pull out the star-forming regions from the underlying 
structure so as to obtain an estimate of the fraction of 
the light of the galaxy that is due to star formation, 
which is basically specific star formation rate. Many fac- 
tors (dust, inclination, resolution, nuclear star formation, 
etc.), however, contribute to create scatter in this simple 
relation, as described in the appendix, so that the cor- 



respondence between clumpiness and sSFR is essentially 
washed out. 

Other morphological parameters are, in principle, also 
affected by the presence of star- forming regions. Asym- 
metry is probably the parameter that is the second most 
sensitive to star formation regions, and it is actually pos- 
sible that the continuous rise of sSFR from A = 0.1 
(where disks start to dominate) to A = 0.4 (where the 
relation becomes flat) ( figure ITB)) be accounted for by the 
presence of star formation regions. On the other hand, 
concentration, Gini and M20, unlike S and A, are also 
sensitive to the presence and size of bulges. The Gini 
coefficient, however, is affected by disk star-forming re- 
gions differently from how concentration and M20 are. 
Compared to a smooth profile, adding bright regions far 
from the center of the galaxy would systematically de- 
crease concentration, but it would, on the contrary, in- 
crease Gini, though star formation by itself can only in- 
crease Gini up to a certain point in the best of cases 
(which is about G ~ 0.5). This behavior means that, 
unlike concentration that can be low even in the pres- 
ence of a bulge due to peripheral bright star forming 
chimps or structures, the Gini coefficient systematically 
increases with the prominence of the bulge. We infer 
that this is the reason for the correlation between the 
Gini coefficient and stellar mass, as bul ges dominate at 
high masses and disks at low masses (|Kauffman n et all 
I2003bt iBundv et al.ll2005l ). The Gini coefficient is, thus, 
among our five parameters, the one that measures struc- 
tural characteristics the most purely, as opposed to mea- 
suring star formation, or a combination thereof. In that 
context, the correlation between star formation rate and 
the Gini coefficient in star-forming galaxies can be sim- 
ply explained by the fact that star formation is, in most 
cases, more intense in nuclear starb ursts and circum nu- 
clear regions than it is in disks (e.g. iKennicu tt 1998). 

Is the fact that different morphological parameters are 
differently affected by young stellar regions a problem for 
our investigation? The answer depends on what we are 
trying to measure with these parameters. We have shown 
in this paper that bulges and disks separate well in planes 
of asymmetry or clumpiness versus a concentration-like 
parameter (C, G and M20), and that that separation is 
fairly independent of the state of star formation in the 
galaxy as, for example, bulges segregate equally whether 
they are old or in the making through intense starbursts. 
Similarly, the asymmetry parameter can be just as large 
for dry mergers as for gas-rich ones. Non-parametric 
morphological parameters thus appear to represent the 
broad, overall structure of galaxies well enough. Their 
usefulness diminishes, however, once one attempts to 
study more detailed morpholgy (e.g. spiral structure, in- 
teraction [other than major merger] signatures, etc.), as 
most parameters can be equally influenced by numerous 
different configurations, thus resulting in scatter in rela- 
tions between morphological and physical parameters, as 
we have discovered throughout this work. 

6. SUMMARY 

In an age of ever expanding ground-based local (SDSS, 
2dF), and spaced-based redshift (COSMOS, GEMS, 
GOODS) surveys, it has become important to be able 
to understand and apply automatic galaxy classification. 
In this paper, we have studied more specifically the con- 
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centration(C), asymmetry (A) and clumpiness(S) param- 
eters as well as of the Gini cocmcicnt(G) and the second 
moment of the brightest 20% of the light (M20) of a sam- 
ple of z ~ 0.7 galaxies in the COSMOS field. 

We observe the galaxy population at redshift z ~ 0.7 
to display a bimodality in both asymmetry and dumpi- 
ness, corresponding to a separation of the bulge-like 
and disk- like galaxies. This partition is further accen- 
tuated when two morphological parameters are taken in 
combination, especially asymmetry or dumpiness with 
a concentration-like parameters (Concentration, Gini or 
M20). Compared to local galaxies, we observe z ~ 0.7 
galaxies to display a similar distribution of morpholo- 
gies, although we see tentative evidence for them to be, 
on average, a little (0 — 50%) more as ymmetric. 

By c ompari ng with local samples o f Kauffmaim et all 
(|2003bft and iBrinchmann et alj (|2004j ). we observe z ~ 
0.7 galaxies to also have similar masses, sizes and sur- 
face mass densities, but higher star formation rates by a 
factor of 0.35 to 0.55 in log SFR. As 0.55 is precisely the 
observed decrease i n star formation rate den sity from 
z ~ 0.7 to present (|Schiminovich et aT1l2005l ). we con- 
clude that this overall decrease in the SFR of disk galax- 
ies, or dimming, accounts for most of the star forma- 
tion density evolution in the Universe in the l ast 6 Gyr. 
We th us corroborate IWolf et ail (|2005[ ) and iBell et afl 
(2005), who reached the same conclusion by respectively 
showing the z ~ 0.7 UV and IR luminosity functions to 
be dominated by normal spiral galaxies, implying that 
the present-day luminosity functions could be mostly 
reached by simple dimming of their studied population. 
Our observed shift in SFR could be, however, on the 
low side (i.e. smaller than 0.55, the overall shift in star 
formation rate density), which suggests some degree of 
number evolution along with a decrease in the individual 
star formation rates. Such a scenario is most consis- 
tent with the observed evol ution of the UV luminosity 
function (|Arnouts et al.ll2005[) . which, besides brighten- 
ing, also includes a moderate steepening of the faint-end 
slope with redshift. 

We also looked at relations between physical and mor- 
phological parameters and find them to possess much 
scatter. We discussed the origin of this scatter for each 
of our morphological parameters, and discover that it is 
multi-faceted. In particular, we find that star formation, 
in certain configurations, can have an effect on derived 
morphological parameters. The resulting combination 
of star-forming regions and underlying structure that is 
measured can thus be, sometimes, hard to interpret. We 
therefore conclude that the set of morphological param- 
eters discussed in this paper, although very useful in de- 
scribing the overall shape of galaxies as demonstrated, 
for example, in the bimodality, is not as well suited to 
study detailed morphology. 

Nevertheless, we find Gini to correlate strongly with 
stellar mass which we understand to be a consequence 
of the fact that the Gini coefficient better traces the 
overall structure of a galaxy than any of the other mor- 
phological parameters we studied. We also observe Gini 
to be the only morphological parameter with which the 
SFR of star- forming galaxies correlates. This correla- 
tion further indicates that many of the strongest star- 
bursts reach values of G similar to that of early-type 
galaxies, which otherwise usually have higher Gini than 



star-forming objects. We see this tendency for starbursts 
(SFR > 10 M Q yr^ 1 ) to carry bulge-like morphologies in 
our other parameters as well, and interpret those obser- 
vations as evidence that episodes of strong star forma- 
tion, proposed as a way to move galaxies from the blu e 
to the red sequence (Bel l et all 12001 iFaber et al.ll2005h . 
are often linked to the growth of a central component. 

We find most of our starbursting galaxies with early- 
type morphologies to have red colors and large UV ex- 
tinction coefficients, implying that they are substantially 
dust-enshrouded. As a consequence, when looking at 
our full sample, we observe a strong correspondence be- 
tween our morphological and FUV - g color bimodalities. 
Because red bulge-dominated galaxies are populated by 
both old and quiescent galaxies on one hand, and dusty 
and star-forming ones on the other, no such strong cor- 
respondence exists between morphology and specific star 
formation rate. Color and morphology are, therefore, 
not necessarily good probes of the physical state of a 
galaxy, but their correlation implies that morphological 
and color evolution predominantly occurs concurrently. 

The most plausible scenario we envisage to explain our 
results is one where bulges grow through episodes of in- 
tense, concentrated, gas and dust-rich episodes of star 
formation, whether from mergers or secular instabilities, 
that confer the galaxies both their red colors and bulge- 
like morphologies as well as their high star formation 
rates. As intense star bursts trigger efficient feedback 
(|Di Matteo et al.l [2005) , the star formation finds itself 
quenched after some time, and the galaxy quickly moves 
to the region of quiescent ellipticals. During that pro- 
cess, both its morphology and color change little since, 
once star formation stops, the FUV-flux rapidly fades (in 
about 300 Myr) and the galaxy remains red. 

Our interpretation implies that a fraction of the stars 
in the red sequence have formed at z < 0.7, and that the 
red sequence has kept growing in the last 6 Gyr through 
episodes of intense star formation, many of which are 
likely linked to merger events. 
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APPENDIX 

EXPLAINING THE NOISY RELATION BETWEEN CLUMPINESS AND SPECIFIC STAR FORMATION RATE 

Intuitively, one would think that dumpiness would be the parameter that best traces the specific star-formation 
rate. However, as we have seen in figures [T51 through 1161 the relation has so much scatter that all correlation, except 
at the low-end (S < 0.12), is washed out. This is confirmed in figure [18] by the histograms of sSFR vs S, as well as 
by the contours of constant specific star formation rate that run largely parallel to the S-axis in all plots. To get a 
better insight into what dumpiness means physically, it is instructive to look at S vs A in figure I17[ which shows 
that sometimes clumpier objects actually have redder colors than less clumpy ones. As mentioned in section 4.1, high 
S/A ratios indicate edge-on galaxies which, because of dust lanes and/or the use of a circular filter in the dumpiness 
algorithm, carry high values of S. They are, however, often symmetric, and thus carry only moderate values of A. Also, 
because of those same dust lanes, they emit less UV light in our direction than would a face-on galaxy. On the other 
hand, some galaxies with blue cores carry diffuse asymmetric features. These galaxies would typically have higher 
asymmetry as well as high UV-fluxes, but would not necessarily have high dumpiness values. These could explain 
why sSFR values are often higher at low S/A ratios and lower at high S/A ratios. They also contribute to reverse 
the expected correlation of sSFR and S. Nevertheless, it is unlikely that inclination effects and diffuse asymmetric 
components are able to entirely account for the wide scatter in the sSFR - S relation. We have seen in section 4.3 that 
dumpiness shows a suspicious correlation with size. It is surprising because all of our morphological parameters are 
normalized quantities. Such a correlation with size can therefore suggest resolution issues, which we now turn to to 
explore. 

Although lElmegreen fc Sa lzcr (1999]) observed that a relatively constant fraction (ss 7%) of the B-band luminosity of 
late spirals and irregulars originates from star-forming complexes, we observe that spiral arms seem to often contribute 
as much to dumpiness as star-forming complexes themselves. As den sity-wave structures are predominantly present 
in large disks, they contribute to create the observed trend. Moreover. lElmegreen et all (|1996f) showed the size of the 
larg est star-forming complexes , that is of the largest clumps, to be more or less proportional to the size of the galaxy, 
and [Elmcgrccn & Salzei] (|1999| ) further showed the number of complexes of a given diameter in a galaxy to be roughly 
proportional to the inverse square of their diameter, while their luminosity goes as the square of their diameter. This 
implies two things. First, we might not be able to resolve any dumpiness in some of our smallest galaxies, and secondly, 
even in galaxies for which we do resolve the largest complexes, since the product n comp i ex x L comp i ex is more or less 
independent of size, the total luminosity we measure in clumps roughly goes as Ljj x {\ogD max — logD m i„}, where 
L d represents a luminosity per dex of size and is multiplied by the difference of the log of the diameter of the largest 
complexes and that of the smallest ones we can resolve. Since our resolution is constant, the total measured luminosity 
of star- forming complexes simply goes as x \ogD max , which in turns is proportional to L ga i axy x log D ga i axy , 
meaning that, since S is normalized to the total flux of the galaxy, S ~ log D ga i axy , or S ~ logr 50 %. The correlation 
between dumpiness and size is thus, at least partly, the consequence of a lack of resolution. This effect adds to the 
flattening of the sSFR y s. S relation as larger disks tend to have lower specific star formation rates (see figure or 
iBrinchmann et al]|2004l ) , and is actually probably an even more significant source of scatter than inclination effects or 
morphological diversity. Because dumpiness is influenced by these many factors, we question its usefulness as a star 
formation indicator. 



